A versatile and interoperable computational framework for the analysis and modeling of COVID-19 disease mechanisms (preprint)

The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing. Community-driven and highly interdisciplinary, the project is collaborative and supports community standards, open access, and the FAIR data principles. The coordination of community work allowed for an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework links key molecules highlighted from broad omics data analysis and computational modeling to dysregulated pathways in a cell-, tissue- or patient-specific manner. We also employ text mining and AI-assisted analysis to identify potential drugs and drug targets and use topological analysis to reveal interesting structural features of the map. The proposed framework is versatile and expandable, offering a significant upgrade in the arsenal used to understand virus-host interactions and other complex pathologies.

Real-world evidence with a retrospective cohort of 15,968 Andalusian COVID-19 hospitalized patients suggests 21 new effective treatments and one drug that increases death risk (preprint)

Despite the extensive vaccination campaigns in many countries, COVID-19 is still a major worldwide health problem because of its associated morbidity and mortality. Therefore, finding efficient treatments as fast as possible is a pressing need. Drug repurposing constitutes a convenient alternative when the need for new drugs in an unexpected medical scenario is urgent, as is the case with COVID-19. Using data from a central registry of electronic health records (the Andalusian Population Health Database, BPS), the effect of prior consumption of drugs for other indications previous to the hospitalization with respect to patient survival was studied on a retrospective cohort of 15,968 individuals, comprising all COVID-19 patients hospitalized in Andalusia between January and November 2020. Covariate-adjusted hazard ratios and analysis of lymphocyte progression curves support a significant association between consumption of 21 different drugs and better patient survival. Contrarily, one drug, furosemide, displayed a significant increase in patient mortality.

Assessing the impact of SARS-CoV-2 lineages and mutations on patient survival (preprint)

After more than two years of COVID-19 pandemic, SARS-CoV-2 still remains a global public health problem. Successive waves of infection have produced new SARS-CoV-2 variants with new mutations whose impact on COVID-19 severity and patient survival is uncertain. A total of 764 SARS-CoV-2 genomes sequenced from COVID-19 patients, hospitalized from 19th February 2020 to 30st April 2021, along with their clinical data, were used for survival analysis. A significant association of B.1.1.7, the alpha lineage, with patient mortality (Log Hazard ratio LHR=0.51, C.I.=[0.14,0.88]) was found upon adjustment by all the covariates known to affect COVID-19 prognosis. Moreover, survival analysis of mutations in the SARS-CoV-2 genome rendered 27 of them significantly associated with higher mortality of patients. Most of these mutations were located in the S, ORF8 and N proteins. This study illustrates how a combination of genomic and clinical data provide solid evidence on the impact of viral lineage on patient survival.

New susceptibility loci for severe COVID-19 by detailed GWAS analysis in European populations (preprint)

Due to the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), deepening the host genetic contribution to severe COVID-19 may further improve our understanding about underlying disease mechanisms. Here, we describe an extended GWAS meta-analysis of 3,260 COVID-19 patients with respiratory failure and 12,483 population controls from Italy, Spain, Norway and Germany, as well as hypothesis-driven targeted analysis of the human leukocyte antigen (HLA) region and chromosome Y haplotypes. We include detailed stratified analyses based on age, sex and disease severity. In addition to already established risk loci, our data identify and replicate two genome-wide significant loci at 17q21.31 and 19q13.33 associated with severe COVID-19 with respiratory failure. These associations implicate a highly pleiotropic ~0.9-Mb 17q21.31 inversion polymorphism, which affects lung function and immune and blood cell counts, and the NAPSA gene, involved in lung surfactant protein production, in COVID-19 pathogenesis.

Real world evidence of calcifediol use and mortality rate of COVID-19 hospitalized in a large cohort of 16,401 Andalusian patients (preprint)

Background: COVID-19 is a major worldwide health problem because of acute respiratory distress syndrome, and mortality. Several lines of evidence have suggested a relationship between the vitamin D endocrine system and severity of COVID-19. Methods: We present a retrospective survival study that includes all Andalusian patients hospitalized between January and November 2020 because of COVID-19 infection. Based on a central registry of electronic health records (the Andalusian Population Health Database, BPS), prescription of vitamin D or its metabolites within 15-30 days before hospitalization were recorded. The effect of treatment with vitamin D metabolites for other indication previous to the hospitalization was studied with respect to patient survival by means of Kaplan-Meyer survival curves and Log Hazard Ratios, using a propensity score to compensate the disbalance of compared classes and the confounding factors. The availability of detailed patient data in the BPS allowed to obtain Real-World Evidence (RWE) of the effects of prior use of vitamin D or its metabolites on the mortality due to COVID-19 infection. Findings: A retrospective cohort of 16.401patients was extracted from the BPS, which includes all the patients hospitalized with COVID-19 diagnosis between January and November 2020 in Andalusia, one of the largest regions in Europe with the size of an average median country. A total of 358 patients were found with cholecalciferol, and 193 with calcifediol, prescriptions 15 days before hospitalization. For a period extended to 30 days before hospitalization, the numbers increase to 416 and 210 and, respectively. Kaplan-Meyer survival curves and hazard ratios support an association between consumption of these metabolites and patient survival. Such association was stronger in calcifediol (Log Hazard Ratio, LHR = -1.27{+/-}0.32) than in cholecalciferol (LHR= -0.56{+/-}0.15), when prescribed 15 days before hospitalization This effect decreases when a larger 30 days period is considered (calcifediol LHR= -1.01{+/-}0.27 and cholecalciferol LHR= -0.27{+/-}0.12), suggesting that the closer was the treatment to the hospitalization the stronger the association. Conclusions: A significant reduction in mortality in patients hospitalized with COVID-19 is associated with the prescription of vitamin D, especially calcifediol, within 15-30 days prior to hospitalization.

From partial to whole genome imputation of SARS-CoV-2 for epidemiological surveillance (preprint)

The current SARS-CoV-2 pandemic has emphasized the utility of viral whole genome sequencing in the surveillance and control of the pathogen. An unprecedented ongoing global initiative is increasingly producing hundreds of thousands of sequences worldwide. However, the complex circumstances in which viruses are sequenced, along with the demand of urgent results, causes a high rate of incomplete and therefore useless, sequences. However, viral sequences evolve in the context of a complex phylogeny and therefore different positions along the genome are in linkage disequilibrium. Therefore, an imputation method would be able to predict missing positions from the available sequencing data. We developed impuSARS, an application that includes Minimac, the most widely used strategy for genomic data imputation and, taking advantage of the enormous amount of SARS-CoV-2 whole genome sequences available, a reference panel containing 239,301 sequences was built. The impuSARS application was tested in a wide range of conditions (continuous fragments, amplicons or sparse individual positions missing) showing great fidelity when reconstructing the original sequences. The impuSARS application is also able to impute whole genomes from commercial kits covering less than 20% of the genome or only from the Spike protein with a precision of 0.96. It also recovers the lineage with a 100% precision for almost all the lineages, even in very poorly covered genomes (< 20%). Imputation can improve the pace of SARS-CoV-2 sequencing production by recovering many incomplete or low-quality sequences that would be otherwise discarded. impuSARS can be incorporated in any primary data processing pipeline for SARS-CoV-2 whole genome sequencing.

Mechanistic modeling of the SARS-CoV-2 disease map (preprint)

Here we present a web interface that implements a comprehensive mechanistic model of the SARS-CoV-2 disease map. In this framework, the detailed activity of the human signaling circuits related to the viral infection, covering from the entry and replication mechanisms to the downstream consequences as inflammation and antigenic response, can be inferred from gene expression experiments. Moreover, the effect of potential interventions, such as knock-downs, or drug effects (currently the system models the effect of more than 8000 DrugBank drugs) can be studied. This freely available tool not only provides an unprecedentedly detailed view of the mechanisms of viral invasion and the consequences in the cell but has also the potential of becoming an invaluable asset in the search for efficient antiviral treatments.

COVID-19 Disease Map, a computational knowledge repository of SARS-CoV-2 virus-host interaction mechanisms (preprint)

We hereby describe a large-scale community effort to build an open-access, interoperable, and computable repository of COVID-19 molecular mechanisms - the COVID-19 Disease Map. We discuss the tools, platforms, and guidelines necessary for the distributed development of its contents by a multi-faceted community of biocurators, domain experts, bioinformaticians, and computational biologists. We highlight the role of relevant databases and text mining approaches in enrichment and validation of the curated mechanisms. We describe the contents of the map and their relevance to the molecular pathophysiology of COVID-19 and the analytical and computational modelling approaches that can be applied to the contents of the COVID-19 Disease Map for mechanistic data interpretation and predictions. We conclude by demonstrating concrete applications of our work through several use cases.

Induced pulmonary comorbidities render CD-1 mice sensitive to SARS-CoV-2 (preprint)

Severe manifestations of COVID-19 are mostly restricted to persons with comorbidities, and they form a significantly high proportion of those which develop life-endangering lung injury. Nevertheless, COVID-19 animal models established to date are not based on preexistence of comorbidities. Here we report that mild pulmonary injury induced by administration of acute-lung-injury stimulants, renders outbred CD-1 mice to be sensitive to SARS-CoV-2. Following intranasal pretreatment of mice with low doses of ricin or bleomycin, SARS-CoV-2 infection caused a severe disease manifested by sustained body loss and mortality rates of >50%. Low-dose-ricin pretreated mice displayed markedly higher levels of viral RNA than mice not pretreated with ricin, not only in the nasal turbinate, trachea and lungs but also in the serum and heart. The deleterious effects of SARS-CoV-2 infection in ricin-pretreated mice were effectively alleviated by passive transfer of polyclonal and monoclonal antibodies generated against SARS-CoV-2 or SARS-CoV-2 RBD. Notably, viral cell entry in the sensitized mice model seems to involve viral RBD binding, albeit by a mechanism other than the canonical ACE2-mediated uptake route. In summary, we present a novel animal model in mice that express native murine ACE2 yet are susceptible to genetically unaltered SARS-CoV-2, for the study of comorbidity-dependent COVID-19 pathology and treatment.

A genome-wide CRISPR/Cas9 knock-out screen identifies the DEAD box RNA helicase DDX42 as a broad antiviral inhibitor (preprint)

Genome-wide CRISPR/Cas9 knock-out genetic screens are powerful approaches to unravel new regulators of viral infections. With the aim of identifying new cellular inhibitors of HIV-1, we have developed a strategy in which we took advantage of the ability of type 1 interferon (IFN) to potently inhibit HIV-1 infection, in order to create a cellular environment hostile to viral replication. This approach led to the identification of the DEAD-box RNA helicase DDX42 as an intrinsic inhibitor of HIV-1. Depletion of endogenous DDX42 using siRNA or CRISPR/Cas9 knock-out increased HIV-1 infection, both in model cell lines and in physiological targets of HIV-1, primary CD4+ T cells and monocyte-derived macrophages (MDMs), and irrespectively of the IFN treatment. Similarly, the overexpression of a dominant-negative mutant of DDX42 positively impacted HIV-1 infection, whereas wild-type DDX42 overexpression potently inhibited HIV-1 infection. The positive impact of endogenous DDX42 depletion on HIV-1 infection was directly correlated to an increase in viral DNA accumulation. Interestingly, proximity ligation assays showed that DDX42, which can be mainly found in the nucleus but is also present in the cytoplasm, was in the close vicinity of HIV-1 Capsid during infection of primary monocyte-derived macrophages. Moreover, we show that DDX42 is also able to substantially decrease infection with other retroviruses and retrotransposition of long interspersed elements-1 (LINE-1). Finally, we reveal that DDX42 potently inhibits other pathogenic viruses, including Chikungunya virus and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2).

Aglycone polyether ionophores as broad-spectrum agents inhibit multiple enveloped viruses including SARS-CoV-2 in vitro and successfully cure JEV infected mice (preprint)

Infections with zoonotic viruses, such as flaviviruses, influenza virus, and the SARS-CoV-2 pandemic coronavirus constitute an increasing global risk. Hence, an urgent need exists for the development of broad-spectrum antivirals to prevent such outbreaks. Here, we show that the maduramycin and CP-80,219 aglycone polyether ionophores exhibit effective broad-spectrum antiviral activity, against various viruses, including Japanese encephalitis virus (JEV), Dengue virus (DENV), Zika virus (ZIKV), and Chikungunya virus (CHIKV), while also exhibiting promising activity against PR8 influenza virus and SARS-CoV-2. Moreover, liposome-encapsulated maduramycin and CP-80,219 provide full protection for mice from infection with JEV in vivo. Mechanistic studies suggest that aglycone polyether ionophores primarily inhibit the viral replication step without blocking endosome acidification to promote the fusion between viral and cellular membranes. The successful application of liposomes containing aglycone polyether ionophores in JEV-infected mice offers hope to the development of broad-spectrum antiviral drugs like penicillin back to 1940s.

Extracellular vesicle-based vaccine platform displaying native viral envelope proteins elicits a robust anti-SARS-CoV-2 response in mice. (preprint)

Extracellular vesicles (EVs) emerge as essential mediators of intercellular communication. DNA vaccines encoding antigens presented on EVs efficiently induce T-cell responses and EV-based vaccines containing the Spike (S) proteins of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) are highly immunogenic in mice. Thus, EVs may serve as vaccine platforms against emerging diseases, going beyond traditional strategies, with the antigen displayed identically to the original protein embedded in the viral membrane and presented as such to the immune system. Compared to their viral and pseudotyped counterparts, EV-based vaccines overcome many safety issues including pre-existing immunity against these vectors. Here, we applied our technology in natural EV's engineering, to express the S proteins of SARS-CoV-2 embedded in the EVs, which mimic the virus with its fully native spikes. Immunizations with a two component CoVEVax vaccine, comprising DNA vector (DNAS-EV) primes, allowing in situ production of Spike harbouring EVs, and a boost using S-EVs produced in mammalian cells, trigger potent neutralizing and cellular responses in mice, in the absence of any adjuvants. CoVEVax would be the prototype of vaccines, where the sole exchange of the envelope proteins on EVs leads to the generation of new vaccine candidates against emerging viruses.

SARS-CoV-2 desensitizes host cells to interferon through inhibition of the JAK-STAT pathway (preprint)

SARS-CoV-2 can infect multiple organs, including lung, intestine, kidney, heart, liver, and brain. The molecular details of how the virus navigates through diverse cellular environments and establishes replication are poorly defined. Here, we performed global proteomic analysis of the virus-host interface in a newly established panel of phenotypically diverse, SARS-CoV-2-infectable human cell lines representing different body organs. This revealed universal inhibition of interferon signaling across cell types following SARS-CoV-2 infection. We performed systematic analyses of the JAK-STAT pathway in a broad range of cellular systems, including immortalized cell lines and primary-like cardiomyocytes, and found that several pathway components were targeted by SARS-CoV-2 leading to cellular desensitization to interferon. These findings indicate that the suppression of interferon signaling is a mechanism widely used by SARS-CoV-2 in diverse tissues to evade antiviral innate immunity, and that targeting the viral mediators of immune evasion may help block virus replication in patients with COVID-19.

In vitro assessment of the virucidal activity of four mouthwashes containing Cetylpyridinium Chloride, ethanol, zinc and a mix of enzyme and proteins against a human coronavirus (preprint)

Background: saliva is established to contain high counts SARS-CoV-2 virus and contact with saliva droplets, contaminated surfaces or airborne particles are sources of viral transmission. The generation of infective aerosols during clinical procedures is of particular concern. Therefore, a fuller understanding of the potential of mouthwash to reduce viral counts and modulate the risk of transmission in medical professional and public context is an important research topic. Method: we determined the virucidal activity of four anti-bacterial mouthwashes against a surrogate for SARS-CoV-2, Human CoV-SARS 229E, using a standard ASTM suspension test, with dilution and contact times applicable to recommended mouthwash use. Results: the mouthwash formulated with 0.07% Cetylpyridinium Chloride exhibited virucidal effects providing a [≥]3.0 log reduction HCoV-229E viral count. Mouthwashes containing 15.7% ethanol, 0.2% zinc sulphate heptahydrate and a mix of enzymes and proteins did not demonstrate substantive virucidal activity in this test. Conclusion: mouthwash containing 0.07% Cetylpyridinium Chloride warrants further laboratory and clinical assessment to determine their potential benefit in reducing the risk of SARS-CoV-2.

Mechanism of SARS-CoV-2 polymerase inhibition by remdesivir (preprint)

Remdesivir is the only FDA-approved drug for the treatment of COVID-19 patients. The active form of remdesivir acts as a nucleoside analogue and inhibits the RNA-dependent RNA polymerase (RdRp) of coronaviruses including SARS-CoV-2. Remdesivir is incorporated by the RdRp into the growing RNA product and allows for addition of three more nucleotides before RNA synthesis stalls. Here we use synthetic RNA chemistry, biochemistry and cryo-electron microscopy to establish the molecular mechanism of remdesivir-induced RdRp stalling. We show that addition of the fourth nucleotide following remdesivir incorporation into the RNA product is impaired by a barrier to further RNA translocation. This translocation barrier causes retention of the RNA 3'-nucleotide in the substrate-binding site of the RdRp and interferes with entry of the next nucleoside triphosphate, thereby stalling RdRp. In the structure of the remdesivir-stalled state, the 3'-nucleotide of the RNA product is matched with the template base, and this may prevent proofreading by the viral 3'-exonuclease that recognizes mismatches. These mechanistic insights should facilitate the quest for improved antivirals that target coronavirus replication.

Ethacridine inhibits SARS-CoV-2 by inactivating viral particles in cellular models (preprint)

More than a million people have now died from COVID-19, because of infection with the SARS-CoV-2 coronavirus. Currently, the FDA has approved remdesivir, an inhibitor of SARS-CoV-2 replication, to treat COVID-19, though very recent data from WHO showed little if any COVID19 protective effect. Here we report that ethacridine, a safe and potent antiseptic use in humans, effectively inhibits SARS-CoV-2, at very low concentrations (EC50 ~ 0.08 M). Ethacridine was identified through a high-throughput screening of an FDA-approved drug library in living cells using a fluorescent assay. Interestingly, the main mode of action of ethacridine is to inactivate virus particles, preventing binding to the host cells. Thus, our work has identified a potent drug with a distinct mode of action against SARS-CoV-2.

Gamma-irradiated SARS-CoV-2 vaccine candidate, OZG-38.61.3, confers protection from SARS-CoV-2 challenge in human ACEII-transgenic mice (preprint)

The SARS-CoV-2 virus caused one of the severest pandemic around the world. The vaccine development for urgent use became more of an issue during the pandemic. An inactivated virus formulated vaccines such as Hepatitis A, inactivated polio, and influenza has been proven to be a reliable approach for immunization for long years. In this pandemic, we produced an inactivated SARS-CoV-2 vaccine candidate by modification of the oldest but the most experienced method that can be produced quickly and tested easily rather than the recombinant vaccines. Here, we optimized an inactivated virus vaccine which includes the gamma irradiation process for the inactivation as an alternative to classical chemical inactivation methods so that there is no extra purification required. Also, we applied the vaccine candidate (OZG-38.61.3) using the intradermal route in mice which decreased the requirement of a higher concentration of inactivated virus for proper immunization unlike most of the classical inactivated vaccine treatments. Thus, the novelty of our vaccine candidate (OZG-38.61.3) is a non-adjuvant added, gamma-irradiated, and intradermally applied inactive viral vaccine. We first determined the efficiency and safety dose (either 1013 or 1014 viral copy per dose) of the OZG-38.61.3 in Balb/c mice. Next, to test the immunogenicity and protective efficacy of the OZG-38.61.3, we immunized human ACE2-encoding transgenic mice and infected them with a dose of infective SARS-CoV-2 virus for the challenge test. We showed that the vaccinated mice showed lowered SARS-CoV-2 viral copy number in oropharyngeal specimens along with humoral and cellular immune responses against the SARS-CoV-2, including the neutralizing antibodies similar to those shown in Balb/c mice without substantial toxicity. This study encouraged us towards a new promising strategy for inactivated vaccine development (OZG-38.61.3) and the Phase 1 clinical trial for the COVID-19 pandemic.

SARS-CoV-2 replication triggers an MDA-5-dependent interferon production which is unable to efficiently control replication (preprint)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third highly pathogenic coronavirus to spill over to humans in less than 20 years, after SARS-CoV-1 in 2002-2003 and Middle East respiratory syndrome (MERS)-CoV in 2012. SARS-CoV-2 is the etiologic agent of coronavirus disease 19 (COVID-19), which ranges from mild respiratory symptoms to severe lung injury and death in the most severe cases. The COVID-19 pandemic is currently a major health issue worldwide. Immune dysregulation characterized by altered innate cytokine responses is thought to contribute to the pathology of COVID-19 patients, which is a testimony of the fundamental role of the innate immune response against SARS-CoV-2. Here, we further characterized the host cell antiviral response against SARS-CoV-2 by using primary human airway epithelia and immortalized model cell lines. We mainly focused on the type I and III interferon (IFN) responses, which lead to the establishment of an antiviral state through the expression of IFN-stimulated genes (ISGs). Our results demonstrate that both primary airway epithelial cells and model cell lines elicit a robust immune response characterized by a strong induction of type I and III IFN through the detection of viral pathogen molecular patterns (PAMPs) by melanoma differentiation associated gene (MDA)-5. However, despite the high levels of type I and III IFNs produced in response to SARS-CoV-2 infection, the IFN response was unable to control viral replication, whereas IFN pre-treatment strongly inhibited viral replication and de novo production of infectious virions. Taken together, these results highlight the complex and ambiguous interplay between viral replication and the timing of IFN responses.

Mechanistic modeling of the SARS-CoV-2 disease map (preprint)

Here we present a web interface that implements a comprehensive mechanistic model of the SARS-CoV-2 disease map in which the detailed activity of the human signaling circuits related to the viral infection and the different antiviral responses, including immune and inflammatory activities, can be inferred from gene expression experiments. Moreover, given to the mechanistic properties of the model, the effect of potential interventions, such as knock-downs, over-expression or drug effects (currently the system models the effect of more than 8000 DrugBank drugs) can be studied in specific conditions. By providing a holistic, systems biology approach to the understanding of the complexities of the viral infection process, this tool will become an important asset in the search for efficient antiviral treatments. The tool is freely available at: http://hipathia.babelomics.org/covid19/